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ABSTRACT 

Cassini-Huygens data are used to re-examine the potential sources of the D /H enhancement over 
solar, measured in methane, in Titan's atmosphere. Assuming that the system is closed with respect 
to carbon, the use of constraints from the Huygens probe for the determination of the current mass 
of atmospheric methane and the most up-to-date determination of D/H from Cassini/CIRS infrared 
spectra allow us to show that photochemical enrichment of deuterium is not sufficient to be the sole 
mechanism yielding the measured D/H value. A possible fractionation between CH3D and CH4 during 
the escape process may slightly enhance the deuterium enrichment, but is not sufficient to explain 
the observed D/H value over the range of escape values proposed in the literature. Hence, alternative 
mechanisms such as a primordial deuterium enrichment must be combined with the photochemical 
enrichment in Titan's atmosphere in order to explain its current D/H value. 
Subject headings: Planets and satellites: formation — Planets and satellites: individual: Titan 



1. INTRODUCTION 

The interpretation of the D/H ratio measured in 
methane in the atmosphere of Titan is not straightfor- 
ward. All measurements indicate values substantially 
higher than the protosolar value, namely the value in 
the hydrogen of the solar nebula, but less than the values 
in water in the Earth's oceans (SMOW) and in comets 
(Bezard et al. 2007). There is no general agreement 
on the cause of the enhancement. Indeed, methane on 
Titan could have originated in the molecular cloud or 
interstellar medium with a high D/H value. Deuterium 
in methane would have then isotopically exchanged with 
molecular hydrogen in the gas phase at the time of en- 
trapment in solids condensed in Saturn's formation zone 
(Mousis et al. 2002) to obtain the D/H value observed 
in Titan. 

On the other hand. Titan could have been accreted 
from solids condensed in an initially hot and dense Sat- 
urn's subnebula. In this picture, the methane incorpo- 
rated in Titan would result from CO and CO2 gas phase 
conversions (Prinn & Fegley 1981) and would present an 
almost solar D/H ratio at the time of its trapping in solids 
ultimately accreted by the satellite (Lunine et al. 1999; 
Mousis et al. 2002). Pinto et al. (1986) and Lunine et 
al. (1999) argued that, because the initial methane reser- 
voir on Titan was likely larger than what is seen today 
in the atmosphere of Titan, the observed D/H enhance- 
ment could be the result of photochemical enrichment of 
deuterium through that isotope's preferential retention 
during methane photolysis (Pinto et al. 1986; Lunine et 
al. 1999). The D/H ratio acquired by the atmospheric 
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methane of Titan would be progressively enriched with 
time via photolysis, until it reaches the value observed 
today. 

Here, we reinvestigate the hypothesis of photochemi- 
cal enrichment of deuterium in the atmosphere of Titan, 
in light of a number of recent Cassini-Huygens measure- 
ments. Pinto et al. (1986) and Lunine et al. (1999) 
estimated the current mass of methane in the atmo- 
sphere of the satellite by assuming that its molar fraction 
is uniform, whatever the altitude. In contrast, we use 
the methane mole fraction and density atmospheric pro- 
files resulting from data collected by the Gas Chromato- 
graph Mass Spectrometer (GCMS) and Huygens Atmo- 
spheric Structure Instrument (HASI) onboard the Huy- 
gens probe during its descent in Titan's atmosphere to 
better constrain the actual mass of methane. Cassini 
data also have provided new constraints on the vertical 
mixing in the atmosphere of Titan, leading Yelle et al. 
(2008) to derive the existence of a CH4 escape flux which 
is about one third the photolytic destruction rate of CH4. 
Here, we investigate the influence of this prodigious es- 
cape on the total fractionation between CH3D and CH4. 

Moreover, we utilize the most recent determina- 
tion of D/H obtained by Bezard et al. (2007) from 
Cassini/CIRS infrared spectra, which indicates values 
of 1.32j^g |^^ X 10"**, substantially higher than those em- 
ployed by Lunine et al. (1999) (7.75±2.25 x 10"^; Orton 
1992). Since the determination of Bezard et al. (2007) 
was obtained by fitting simultaneously the vq bands of 
both 13CH3D and i^CHgD and the U4 band of CH4 from 
precise information on the CH4 mixing ratio and tem- 
perature profile in the stratosphere available from the 
Huygens descent and limb-viewing Composite Infrared 
Spectrometer (CIRS) measurements, we believe that this 
measurement is more reliable than the previous ones. 

All these revisions, together with the use of updated 
rate coefficients for methane loss derived from a recent 
photochemical model (Vuitton et al. 2008), allow us to 
show that the photochemical enrichment of deuterium is 
not efficient enough in the atmosphere of Titan to ex- 
plain its current D/H value, even if the current atmo- 
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spheric reservoir of methane is postulated to exist over 
4.5 Gyr. We conclude that the D/H ratio in methane 
initially trapped in Titan was already higher than the 
protosolar value prior to its release in the atmosphere. 

2. ISOTOPIC ENRICHMENT MODEL 

Our photochemical enrichment model is derived from 
the formalisms of Pinto et al. (1986) and Lunine et al. 
(1999). Wc determine the deuterium enrichment that oc- 
curred following the ultimate outgassing event that gave 
existence to the current atmospheric methane of Titan. 
This event might have started 4.5 Gyr ago at the epoch 
of core overturn, giving birth to a reservoir of methane 
dense enough to allow its survival up until now. On the 
other hand, the current reservoir of methane might also 
have resulted from a much more recent outgassing event, 
as late as ~ 0.6 Gyr ago based on evolution models (To- 
bie et al. 2006) which are consistent with the surface 
cratering record (Lorenz et al. 2007). We assume that 
the system (atmosphere and subsurface from which the 
methane outgasses) is open with respect to the dissoci- 
ated hydrogen (either deuterium or protium) which is 
lost to space, and to carbon only when methane loss is 
considered in the upper atmosphere of Titan. 

We define R as the ratio of the total mass of methane 
expelled from the interior of Titan (and constituting 
the initial reservoir) to the current atmospheric mass of 
methane. Let be Ni{t) and N2{t) the respective total 
column abundance of CH4 and CH3D on Titan at time 
t. With the condition A^i > iV2, the evolution of Ni{t) 
and N2{t) is given by 



dNi 
"dT 



diVz 
dt 



Niit) 



m{t)+N2it) 
-(F + $) 



F 



mit) + N2it) 



(1) 



N2{t) 



^N,{t) 
'{qF + l<i>) 



N2{t) 
N2{t) 

Ni{t) 



F-l- 



N2{t) 



Niit)+N2it) 



(2) 



where F is the net photolytic destruction rate of CH4 and 
$ the methane escape rate. The parameter q corresponds 
to the ratio of k2 over ki, namely the respective rates 
for CH3D and CH4 destruction, and ranges between 0.8 
and 0.88 from consideration of the chemical kinetics of 
deuterated species (Lunine et al. 1999). The value of the 
free parameter I ranges between and 1, representing a 
possible fractionation between CH3D and CH4 during 
the escape process. 

The integration of Eq. [T]and Eq. [2] leads to the follow- 
ing expression for R: 
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where / is defined as the ratio of D/H observed in Ti- 
tan's current atmospheric methane to protosolar D/H 
((D/H)0 = 2.35 ± 0.3 X 10"^ Mousis et al. 2002). 

Alternatively, R can be expressed as follows (Lunine et 
al. 1999): 
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where g is the Titan's surface gravity, rncH4 the mass of a 

( s) 

methane molecule, 2:^4 the methane surface mole frac- 
tion (a;[f^^ = 4.9 x 10"^; Niemann et al. 2005) and P'-"^ 
the pressure at the surface of Titan. The time elapsed 
since the formation of the initial methane reservoir up to 
now is T. Since the term x^j^ P*^"-* /g corresponds to the 
total mass of methane per unit of area, R can be written 
as: 
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where Mcn,^ is the cumulative mass of atmospheric 
methane per unit area. 

We adopt here F = 6.9 x 10^^ molecules m~^.s~^, a 
value derived from the photochemical model of Vuitton 
et al. (2008). In addition, when the escape of CH4 is 
considered in our calculations, we set $ = 2.8 x 10^'^ 
m^^.s^^, namely the value of CH4 escape derived by 
Yelle et al. (2008) from Cassini INMS measurements. 
Otherwise, $ is set to and the deuterium enrichment 
is calculated in a way similar to the work of Pinto et al. 
(1986) and Lunine et al. (1999). 

Huygens probe data allow us to accurately determine 
the value of Mch,- Indeed, knowing the atmospheric 
density and methane mole fraction profiles (see Figs. 1(a) 
and 1(b)), one can write: 
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where Xch, is the methane mass fraction, p(z) the den- 
sity profile and z the altitude. H is the altitude above 
which Men, has no significant change. Figure[T]^c) shows 
that McH4 — 2112 kg.m^^ is asymptotically reached 
when the altitude exceeds ^ 70 km. This determination 
is subtantially lower than the value of ~ 5338 kg.m~^ 
derived from Lunine et al. (1999). This discrepancy re- 
sults from the use by Lunine et al. (1999) (as well as by 
Pinto et al. 1986) of a uniform methane mole fraction 
in the atmosphere of Titan whatever the altitude, in the 
absence of available in situ measurements. 

3. RESULTS 

We consider two different cases of photochemical en- 
richment. In the first case, the CH4 escape from the 
upper atmosphere of Titan is neglected. In the second, 
this escape and a consequent additional fractionation be- 
tween CH3D and CH4 are included. 

Figure [2] summarizes the results for deuterium enrich- 
ment via photodissociation calculated with Eq. [3] assum- 
ing that $ = 0, and shows the initial methane reservoir 
R (normalized to the present one) against q. Three cases 
of present-day deuterium enrichment in the atmospheric 
methane of Titan are represented: the solid curve cor- 
responds to the nominal value reported by Bezard et 
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TABLE 1 

Deuterium enrichment expected in the methane of Titan's 

CURRENT atmosphere. 
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Fig. 1. — (a) Density p of Titan's atmosphere as a function 
of the altitude z derived from the HASI data (Fulchignoni et al. 
2005). (b) Methane mole fraction xch^^ as a function of the altitude 
derived from the GCMS data (Niemann et al. 2005). Diamonds 
show the in situ measurements and the solid line corresponds to the 
values interpolated with the spline technique, (c) Cumulative mass 
of methane in Titan's atmosphere as a function of the altitude. 

al. (2007) (/ — 5.6) and the dashed curves to the ex- 
treme vahies (/ = 4.5-7.2) obtained when uncertainties 
are taken into account. We assume in all our calculations 
that the D/H ratio in the methane initially acquired by 
Titan during its accretion is protosolar. 

Note that a smaller value of q yields greater fractiona- 
tion for a given amount of methane photolysis, because 
it corresponds to deuterium being more tightly bound 
(Lunine et al. 1999). Therefore, for a given deuterium 
enrichment, smaller values of q allow the initial methane 
reservoir of Titan to be smaller than in the case required 
for higher values of q. The two horizontal lines repre- 
sent values that would be acquired by R if the actual 
methane reservoir were to exist since 0.6 or 4.5 Gyr ago, 
respectively (see Eq. [S]). Figure [2] shows that the ini- 
tial reservoir of Titan's atmospheric methane was ~ 16 
times more massive than the current one if it was formed 
0.6 Gyr ago, provided that F remained fixed to its cur- 
rent value throughout the existence of this reservoir. If 
methane were present in the atmosphere of Titan since 
4.5 Gyr ago, the initial reservoir was ~ 126 times more 
massive than the current one. 

Table [T] summarizes the deuterium enrichments via 
photolysis in the atmospheric methane calculated for the 
adopted limits on plausible values of q and the two dif- 
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Fig. 2. — The fractionation of deuterium in methane photochem- 
istry with $ = 0. Top left: plotted is ij as a function of q (see Eq. 
[Sjl. Three curves are shown, corresponding to different present-day 
deuterium enrichments measured in methane by Cassini/CIRS (see 
text). Bottom: the two horizontal lines represent values that would 
be acquired by R for two different values of r (see Eq. |5)l . The two 
vertical lines represent limits on plausible values of q. 

ferent values of r, in the case $ = 0. This table shows 
that, assuming a protosolar D/H in the methane orig- 
inally released into the surface-atmosphere system, the 
photochemical enrichment is not efficient enough to al- 
low the atmospheric D/H to reach the observed enrich- 
ment, even if the reservoir is postulated to have existed 
since 4.5 Gyr. A higher D/H ratio than the protosolar 
value must be advocated in the methane of Titan prior 
its outgassing, in order to explain the observed enrich- 
ment. Depending on the adopted value for g, the range 
of initial deuterium enrichments /o needed by the initial 
methane reservoir to allow photolysis to reach the nom- 
inal value of Bczard et al. (2007) is between 3.2 and 4.0 
for T — 0.6 Gyr, and between 2.2 and 3.2 for r — 4.5 
Gyr. 

Figure [3] shows the results for deuterium enrichment 
via photodissociation calculated with Eq. [3] assuming 
that <i> = 2.8 X 10^^ m~^.s~^, and shows R against q. The 
two upper left curves have been determined for / = 5.6 in 
Eq. [31 namely the nominal value reported by Bezard et 
al. (2007). The solid curve corresponds to the case where 
no fractionation occurs between GH3D and CH4 during 
methane escape, while the dot-dot-dashed curve assumes 
there is fractionation. A word of caution must be given 
here. We have arbitrarily set the escape rate equal to the 
CH4/CH3D molar ratio (/ — 16/17) in order to quantify 
the influence of a possible fractionation between these 
two species. This is an order-of-magnitude estimate be- 
cause the escape mechanism of methane remains unclear 
(Yelle et al. 2008). As a result, the fractionation may be- 
have differently. Figure [3] shows that, despite the higher 
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Fig. 3.— Same as Fig. [2] but with <!> = 2.8 X 10^^ m-^.s-l. Top 
left: both curves show ij as a function of q in the case of / = 5.6. 
The sohd curve corresponds to the case where no fractionation is 
assumed during the methane escape (Z = 1). The dot-dot-dashed 
curve corresponds to the case where a fractionation between CH3D 
and CH4 occurs. 

values of R obtained with the methane loss in the upper 
atmosphere of Titan {R ~ 22 for r = 0.6 Gyr and R ~ 
166 for T — 4.5 Gyr), the discrepancy between the values 
of R{q) and R{t) increases when compared with Fig. [21 
even if a fractionation between CH3D and CH4 is taken 
into account. This translates into a lower photochemical 
enrichment than when methane escape is neglected (see 
Table [1]). The range of /o needed by the initial methane 
reservoir to allow the matching of photolysis with the 
nominal value of Bezard et al. (2007) is consequently 
higher than without methane escape. Here, /o is be- 
tween 3.6 and 4.3 for t = 0.6 Gyr, and between 2.7 and 
3.6 for r = 4.5 Gyr, depending on the adopted value of q 
and with no fractionation. If the fractionation between 
CH3D and CH4 is considered, the range of /o values is 
slightly lower, but still higher than the one inferred in 
the case of no methane escape. 

4. DISCUSSION 

Our results differ from those obtained by Lunine et 
al. (1999) who concluded that the deuterium enrich- 
ment via photolysis was efficient enough to explain the 
current D/H value observed in Titan's atmosphere even if 
the methane incorporated in the forming satellite had ac- 
quired only a slightly supersolar D/H. By considering the 
Cassini-Huygens data, we show that the minimum value 
required for /o is higher than in Lunine et al. (1999), 
even when the escape of methane and a possible fraction- 
ation between CH3D and CH4 are included. Therefore, 
additional mechanisms must be considered in order to 
explain the D/H value observed in Titan's methane. 



One possibility is isotopic thermal exchange of molec- 
ular hydrogen with CH3D in the gas phase of the solar 
nebula. In this case, CH3D might have originated from 
highly deuterium-enriched ices that vaporized when en- 
tering the nebula (Mousis et al. 2002). During the cool- 
ing of the nebula, methane was trapped in crystalline 
ices around 10 AU, perhaps clathrates formed at ^ 60 
K, and incorporated into planetesimals that were pre- 
served in the Saturnian subnebula because the latter was 
too cold to fully vaporize the ices (Mousis et al. 2002; 
Ahbert & Mousis 2007). Subsequent to outgassing, ad- 
ditional enrichment in D/H would have occurred by the 
photochemical and escape mechanisms quantified here. 

Atreya et al. (2006) have argued that methane formed 
in Titan's interior by reaction of CO and CO2 with water 
and rock, in which case the D/H ratio in methane out- 
gassed to the atmosphere would have been determined 
by that in the water, modified by fractionation associated 
with these so-called "serpentinization" reactions. Our re- 
sults show that such a model for methane's origin must 
be able to produce a substantially supersolar value of 
D /H in the outgassed methane. Computation of the D /H 
fractionation associated with serpentinization depends 
on the details of the temperature in Titan's interior, the 
original D/H in the water, and possibly on the compo- 
sition of the silicate buffer. It is unlikely to be large in 
view of the requisite high temperatures well above the 
water melting point) in Titan's interior. Thus, for this 
case, the considerations we made above for primordial 
methane would have to be invoked instead for the wa- 
ter, namely, that it was significantly enriched in D/H in 
the planetesimals that formed Saturn's largest moon. Fi- 
nally, because the original appearance of methane in Ti- 
tan's atmosphere could have been quite early, some 4.5 
billion years ago, its origin might be tightly connected to 
the formation of the N2 atmosphere from the dissociation 
of primordial NH3. Hence, important amounts of H and 
H2 could be present in the ancient atmosphere, leading to 
the possibility of additional fractionation between CH4 
and CH3D via the following reversible equations CH4 -|- 
HD ^ CH3D + H2 and C + CH3 = H + CH2D (see Lee 
et al. 2001 for a list of reactions forming back CH4 and 
CH3D from CH3). However, it is difficult to quantify the 
influence of these reactions on the final CH3D/CH4 ratio 
because the existence of a H2-rich primordial atmosphere 
remains uncertain. 
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sity profile of Titan as measured by the Huygens probe. 
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